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We propose to apply Ramsey's method of separated oscillating fields to the spectroscopy of the 
quantum states in the gravity potential above a vertical mirror. This method allows a precise mea- 
surement of quantum mechanical phaseshifts of a Schrodinger wave packet bouncing off a hard 
surface in the gravitational field of the earth. Measurements with ultra-cold neutrons will offer a 
sensitivity to Newton's law or hypothetical short-ranged interactions, which is about 21 orders of 
magnitude below the energy scale of electromagnetism. 

PACS numbers: 03.65.Ge,03.65.Ta,04.50.-h,04.80.Cc,11.10.Kk 
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INTRODUCTION 



The system of a Schrodinger quantum particle with 
mass m bouncing in a linear gravitational field is known 
as the quantum bouncer [J, 0, Q and "Quantum wave 
packet revivals" can be found in Gravity tests with 
neutrons as quantum objects or within the classical limit 
are reviewed in Above a vertical mirror, the linear 
gravity potential leads to discrete energy eigenstates of 
a bouncing quantum particle. Such quantum states have 
been demonstrated at the Institut Laue-Langevin with 
ultra-cold neutrons in a previous collaboration [f| 0, HI ■ 
Above a vertical mirror, linear gravity potential leads to 
discrete energy eigenstates of a bouncing particle. The 
lowest energy eigenvalues E n , (n = 1, 2, 3, 4, 5), are 1.41 
peV, 2.46 peV, 3.32 peV, 4.09 peV, and 4.78 peV. The 
energy levels together with the neutron density distribu- 
tion are shown in Fig. [TJ The idea of observing quantum 
effects in such a gravitational cavity was discussed with 
neutrons [l(j or atoms 

An important feature of the quantum bouncing ball - 
in contrast to the harmonic oscillator problem - is the fact 
that levels are not equidistant in energy. A combination of 
any two states can therefore be treated as a two-level sys- 
tem. The energy eigenstates in the gravity potential can 
be coupled to a mechanical or magnetic oscillator field. 
Transitions between quantum states in the gravitational 
field of the earth i.e. a change of the state occupation 
can therefore be induced similar to magnetic transitions, 
which occur when the oscillator frequency equals one of 
the Bohr frequencies of the system. This magnetic reso- 
nance method was in the original conception for measure- 
ments of nuclear magnetic moments [HI, [l3| , but soon it 
became a very general technique for radiofrequency spec- 
troscopy [lij . Ramsey developed his method of separated 
oscillating fields in which the oscillatory field is confined 
to a region at the beginning and a region at the end with 
no oscillating field in between [l5| . Variations of Ramseys 
method is inherently connected with precision measure- 
ments ranging from atomic clocks [l6[ to atom interfero- 
metry [l7[, from NMR [l8| to quantum- metrology [lj|, 
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FIG. 1: Energy eigenvalues (blue) and neutron density distri- 
butions (red) for level one to five. 



or the related spin-echo technique [20(. That method has 
also been used to measure the precession frequency of 
atoms, molecules or neutrons in a weak magnetic field, 
for example in a search for permanent atomic or neu- 
tron electric-dipole moments and in constructions of sen- 
sitive magnetometers. The sensitivity is extremely high, 
because a quantum mechanical phase shift is converted 
into a frequency measurement. The sensitivity reached so 
far [2l[ in a search for the electric dipole moment of the 
neutron is 6.8 x 10 -22 eV, or one Bohr rotation within 6 
days. 

In analogy to these examples from electrodynamics, 
we discuss here an application of Ramsey's method to 
probe the eigenstates in the gravity potential. Such a tech- 
nique should open a new way to precision gravity experi- 
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ments and we propose to apply it to quantum states of 
neutrons or atoms in the gravitational field of the earth. 
This method will allow a precise measurements of energy 
differences with a precision similar to the magnetic res- 
onance technique. Here, we are sensitive to energy shifts 
of a Schrodinger wave packet bouncing off a hard sur- 
face. Such energy shifts are expected from hypothetical 
gravity-like forces in the light of recent theoretical devel- 
opments in higher-dimensional field theory and will al- 
low searches for pseudo-scalar coupling of axions in the 
previously experimentally unaccessible astrophysical ax- 
ion window 122,531, see Sec. III. 



II. RAMSEY'S METHOD AND ITS 
APPLICATION TO GRAVITY POTENTIALS 



A quantum mechanical system that is described by 
two states can be understood in analogy to a spin 1/2 sys- 
tem (assuming two states of a fictitious spin in the multi- 
plet, similarly to spin up and spin down states). The time 
development of such systems is described by the Bloch 
equations. In magnetic resonance of a standard spin 1/2 
system, the energy splitting results in the precession of 
the related magnetic moment in the magnetic field. Tran- 
sitions between the two states are driven by a transverse 
magnetic radio frequency field. Similar concepts can be 
applied to any driven two level system, e.g. in optical 
transitions with light fields. Here we apply this picture to 
quantum states in the gravity field. 

We start with a short description of Rabi's 
method [12| to measure the energy difference between a 
two-level system with a coupled oscillating field. With 
LOp q , the frequency difference between the two states, ui, 
the frequency of the driving field, f2jj, the Rabi frequency 
and the time t, the Hamiltonian H is given by 



As we will show below, one can drive transitions be- 
tween quantum states in gravity above the mirror by vi- 
brating the mirror surface. 

Lets consider the motion of ultracold neutrons in the 
gravitational field above a mirror. We assume the gravi- 
tational force to act in -z-direction, while the mirror is 
aligned with the xy-plane, vibrating with amplitude a 
in z-direction. The motion in x- and y-direction is free 
and completely decouples from that in z-direction. It suf- 
fices therefore to consider the time-dependent Schrodinger 
equation restricted to the z-direction 



2m dz 2 



mgz + VqQ(— z + asmcut) > & = ih 



<9* 
~dt 



(4) 

Here, g is the acceleration of gravity, m is the mass of the 
neutron and 6 is the Heaviside step function. The poten- 
tial Vb « 100 neV associated with the substance of the 
mirror is repulsive and much larger than eigenenergies 
of the lowest quantum states in the gravitational field. 
Therefore Eq. (0| must be solved with the boundary con- 
dition ^(z = asmu>t,t) = 0. For further considerations it 
is preferable to introduce z = z — asmuit and to transform 
Eq. into the rest frame of the mirror, 



where 



9* 

{H + W(z,t)}* = to- 
ot 



h 2 d 2 

H = - — -xzz + mgz + Vq®(-z) . 
2m oz 1 



W(z,t) = a 
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and &(z, t) — &(z, t). The hamiltonian Hq describes the 
neutron in the gravitational field above a mirror at rest. 
The second term W(z, t) accounts for the vibration of the 
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The solution can be expressed in terms of the eigen- 
(1) functions ip n (z) of H 



The probability of being found in the excited state 
as a function of time is 

p (" = (l) 2sta2 (^)' (2) 

where the effective Rabi frequency is 

n- R = ^Jni + iujpg-Lu) 2 = ^Jn 2 R + 6 2 , (3) 

with detuning S from resonance. The sinusoidal popula- 
tion transfer is referred to as Rabi flopping. It has been 
proposed to measure the energy levels of a neutron in 
the gravitational field of the Earth with this method 
(GRANIT experiment [13, HE] ). The periodic drive is given 
by neutrons moving through a spatially oscillating mag- 
netic field created by horizontal conducting wires. 



*(z,t) = ^C„(t)e- lB " t /VnW (8) 

n 

with time-dependent coefficients C n (t). Projection of 
Eq. ([5]) on the eigenstates of Hq yields a system of differ- 
ential equations for the coefficients C n (t) 

ih ■ jC n (t) = Ys^nlW^k) ■ C k (t) ■ . (9) 

k 

The transitions between different quantum states is gov- 
erned by the matrix elements of W(z, t) defined in ([7]) 

{tp n \W\ipk) = a [mgS n> k sin u>t + ihu>Q nt k cos cot] (10) 

with 

f°° d 
Qn,k = / dzip n (z)— i> k (z) . (11) 
J dz 
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TABLE I: Relevant overlap integrals Q n ,k defined in Eq. (|11[) 
for the five lowest eigenstates in the gravitational field in 



The relevant overlap integrals Q n ,k for the transitions be- 
tween the lowest eigenstates in the gravitational field are 
given in Tab. |TJ 

The physics behind the transitions between the en- 
ergy eigenstates of the quantum bouncer caused by a vi- 
brating mirror or an oscillating potential is related to ear- 
lier studies of energy transfer when matter waves bounce 
of a vibrating mirror |26l, |27l l28l . |29| or on a time depen- 
dent crystal [30l . l3ll . |32|. In the later cases the transitions 
are between continuum states, in the quantum bouncer 
between discrete eigenstates. Most interesting for our pro- 
posal to drive transitions between eigenstates of the quan- 
tum bouncer with a vibrating mirror is the physics of 
reflection of a neutron by an oscillating potential step as 
has been investigated at the research reactors Munich and 
Geesthacht (FRM and FRG) [HI, however in a different 
energy regime. 

Applying Ramsey's resonance method with sepa- 
rated oscillating fields will allow a careful measurement of 
the energy eigenstates states of the quantum bouncer [33| . 
We propose to implement it with neutrons by traversing 
five regions as shown in Fig. [5] The horizonal direction 
in space is considered as free motion, while the verti- 
cal one is described by a one-dimensional time dependent 
Schrodinger equation (see e.g. Eq. |4|). 

To implement Ramsey's method, one has to realize 
(1) a state selector, (2) a region, where one applies a it/2 
pulse creating the superposition of the two states, whose 
energy difference should be measured, (3) a region, where 
the phase evolves, (4) a second region to read the relative 
phase by applying a second ir/2 pulse, and finally (5) a 
state detector. 

In the following we will describe all these components 
as they are shown in Fig. [21 

In region one, neutrons are prepared in a specific 
quantum state \p) in the gravity potential following the 
procedure demonstrated in Q . A polished mirror on bot- 
tom and a rough absorbing scatterer on top at a height 
of about 20 /jm is a realization of a state selector. It pre- 
pares neutrons into the ground state. Neutrons in higher, 
unwanted states are scattered out of the system and ab- 
sorbed i.e. C\ = 1 and C n = for n > 1. A quantum 
mechanical description of such a system can be found 
in Q. The neutron passage through a mirror-scatterer 
system has also been studied in a frame, where the rough 
scatterer surface has been treated as a time-dependent 




FIG. 2: Sketch of the proposal. Region 1: Preparation in a 
specific quantum state, e.g. state one with polarizer. Region 2: 
Application of first 7r/2-flip. Region 3: Flight path with length 
L. Region 4: Application of second 7r/2-flip. Region 5: State 
analyzer. 



variation of the scatterer position [34( . 

In region two of length I, the first of two identical 
oscillators is installed. Here, transitions between quantum 
states \p) and \q) are induced within a time r according 
to Equ. [2j The oscillator frequency at resonance for a 
transition between states with energies E q and E p is 

_ (Eg - E p ) 

Upq~ ^ • 

The squared ratio of and JTr as a function of the 
driving field u for transitions |1) — > |3), |1) — > |2) and |2) 
— » |3) is shown in Fig. [3J For the |1) — * |2) transition, 
which we chose as an example for transitions in a two 
level system, u>i2 = u>2 — u>i = 2tt x 254 s _1 . 
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FIG. 3: Ratio as a function of lj [s _1 ] for transition 

|1 >^ |3 >, |1 >^ [2 > and [2 >^ |3 >. The strength of the 
vibration was set to 1.0 m/s 2 . 
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Driven on resonance (to 



j), this oscillator drives 



the system into a coherent superposition of state \p) and 
\q). A 7r/2-pulse, that is one with pulse area Qrt = n/2, 
creates an equal superposition between state \p) and \q). 
This can be done by using oscillating magnetic gradient 
fields or by vibrating mirrors i.e. a modulation of the mir- 
ror potential in height. 

In the intermediate region three, a non-oscillating 
mirror with a neutron flight path of L and flight time T 
follows. It might be convenient to place a second mirror 
on top of the bottom mirror at a certain height h. It allows 
to tune the resonance frequency between \a) and \b) due 
to the additional potential, and it provides an effective 
doubling of sensitivity in a search for hypothetical axion 
induced phase shifts or other fifth forces, see section HTT1 

Subsequently, in region four a second oscillator in 
phase with the oscillator in region two is placed. If the 
oscillating u) is equal uj pq than the system is at resonance 
and we have a complete reversal of the state occupation 
between \p) and \q). There is no change in the relative 
phase of the oscillator and the quantum state of the neu- 
tron independent of the neutron velocity. In the other 
cases for u> ^ aj pq , a velocity dependent relative phase 
shift builds up, since a slower neutron is in the region 
three longer and experiences a greater shift than a faster 
neutron. 

Afterwards in section five, such a phase shift can be 
measured by transmission through a second state selector. 

This method can be realized with some modifications 
to the previous setup in the following way: Neutrons are 
taken from the ultracold neutron installation PF2 at ILL 
with a measured horizontal velocity n = 3.2 m/s < « < 20 
m/s. At the entrance of the experiment, a collimator ab- 
sorber system limits the transversal velocity to an energy 
in the pico-eV range. The experiment itself is mounted on 
a polished plane granite stone with an active and passive 
antivibration table underneath. This stone is leveled us- 
ing piezo translators 35]. Inclinometers together with the 
piezo translators in a closed loop circuit guarantee level- 
ing with a precision better than 1 /irad [36| . A solid block 
with dimensions 10 cm x 15 cm x 3 cm composed of op- 
tical glass serves as a mirror for neutron reflection. The 
neutrons see a surface that is essentially flat. In region 
one, an absorber/scatterer that is a rough mirror with 
a surface roughness of about 0.4 ^m is placed above the 
first mirror at a height of 27 pm in order to select the first 
quantum state. The other states are efficiently removed, 
except for the second state, which is still present with a 
contribution of a few percent. In region two, a second mir- 
ror is placed after the first one. Piezo elements attached 
underneath induce a fast modulation of the surface height 
with amplitude a according to Eq. 2] 

As an example, we consider transitions between state 
|1) -> |2) for the most probable velocity at the PF2/UCN 
beam position, 6 m/s. The length I = 15 cm of this mir- 
ror is chosen in such a way to provide a neutron in a 
superposition of these two quantum states after r = 25.0 



ms. Region three has a flight path of L = 80 cm on a 
single mirror between the two oscillators in region two 
and, identical to region two, in region four. In region five, 
a state selector as an analyzer is placed, identical to the 
selector in region one but with a neutron detector behind 
for counting the transmitted neutrons. Calculated tran- 
sition probabilities for \p) and \q) as a function of 
lo is shown in Fig. [4] for different parts of the measured 
velocity spectrum. 

This method can also be applied to stored ultra-cold 
neutrons. Fig. [5] shows the theoretical Ramsey signal for 
a neutron storage time of 100 s. The appeal of a neutron 
storage lies in a very narrow resonance line. A search for 
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FIG. 4: Transition probability for a neutron velocity 5.5 m/s 
< v < 6.5 m/s (red) and 3.0 m/s < v < 9.0 m/s (blue) as a 
function of uj [s -1 ]. The strength of the vibration was set to 
1.0m/s 2 . 
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FIG. 5: Theoretical Ramsey signal for a neutron storage time 
of 100 s. Again, the strength of the vibration was set to 
1.0m/s 2 . 



phase shifts are suggested in the next section. 
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III. PHASE SHIFTS FROM HYPOTHETICAL 
GRAVITY-LIKE FIFTH FORCES 

Theoretical considerations arising from higher- 
dimensional gravity, gauge forces or massive scalar fields 
suggest that the Newtonian gravitational potential for 
masses irii and rrij and distance r should be replaced by 
a more general expression including a Yukawa term, 

V( r )=_G^^(l-a.e- r / A ), (13) 
r 

where A is the Yukawa distance over which the corre- 
sponding force acts and a is a strength factor in units of 
Newtonian gravity. G is the gravitational constant. Most 
interesting, from the experimental point of view, are sce- 
narios, where the strength of the new force is expected to 
be many orders of magnitude stronger than Newtonian 
gravitation. Such forces a re p ossible via abelian gauge 
fields in the bulk [S3, 0, M, E3 , (see also [U El for 
explicit realizations in string theory) . The strength of the 
new force would be 10 6 < a < 10 12 stronger than grav- 
ity, independent of the number of extra dimensions n [3^ | . 
The observation of quantum states already tests specula- 
tions of this kind on large extra dimensions of submillimc- 
ter size of space-time [22|, 0, |44| . Most recent theoretical 
developments support the original proposal of large ex- 
tra dimensions with bulk gauge fields and more specific 
predictions for a high interaction strength can be made. 
One proposal of Callin et al. [45| predicts deviations from 
Newton's law on the micron scale on the basis of super- 
symmetric large extra dimensions (SLED). The basic idea 
behind this proposal is to modify gravity at small dis- 
tances in such a way as to explain the smallness of the 
observed cosmological constant. The hope is to achieve 
this without changing non-gravitational physics (which 
may be possible because of the small gravitational re- 
sponse of the vacuum in specific models of SLED) and 
to link the size of the extra dimensions to the energy 
density 10 10 eV/m 3 governing the observed dark energy 
component of the universe. In the concrete constructions 
presently under discussion a radius r of 10 microns as well 
as the necessary interaction strength up to a = 10 6 may 
turn out to be well-motivated. 

Other experimental limits on extra forces are de- 
rived from mechanical exp eriments and can be found, e.g., 
in [H, H3, IH, IM H3, Hit HI- Probing sub- micron forces 
by interferometry of Bosc-Einstein condensated atoms has 
been proposed by [54j . In practice, the experimental data 
are subject to corrections, which can be orders of mag- 
nitude larger than the effects actually searched for. It is 
therefore important to stress the completely different na- 
ture of possible systematic effects inherent to these micro- 
mechanical experiments as compared to those in neutron 
experiments. In the former case gravitational interactions 
are studied in the presence of large van der Waals and 
Casimir-forces, which depend strongly on the geometry 
of the experiment, and the theoretical treatment of the 
Casimir effect is a difficult task. Currently, atomic force 
microscopes (AFM) measurements using functionalized 



tips determine the limits on non-Newtonian gravitation 
below 10 pm. The best experimental data available ob- 
tained with AFM, are claimed to be at the same level 
of accuracy (1-2%) as the numerical calculations of the 
Casimir force. The major obstacle for improvement in the 
theoretical calculation is the fact that it is very hard to 
take the boundary conditions of the tip and its function- 
alization properly into account. 

Our approach of probing Newtonian Gravity at the 
micron scale with the help of Ramsey's Method of Separa- 
ted Oscillating Fields is advantageous because of its small 
systematic effects. In contrast to atoms the electrical po- 
larizability of neutrons [53| inducing such Casimir effects 
or van der Waals forces is extremely low. This together 
with its electric neutrality the neutron provides the key to 
a sensitivity of more than 10 orders of magnitude below 
the background strength of atoms. 

The dynamics of such a quantum mechanical wave 
packet combines quantum theory with aspects of Newto- 
nian mechanics at short distances. When a neutron with 
mass m approaches the mirror, the mass of this extended 
source might modify the earth acceleration g, when strong 
non-Newtonian forces with range A and strength a are 
present. For small neutron distances z from the mirror, 
say several micrometres, we consider the mirror as an in- 
finite half-space with mass density p. By replacing the 
source mass mi by drrii and integrating over drrii, the 
modified Newtonian potential AV(z) is having the form 

AV(z) = 2irmpa\ 2 Ge- lzl/x 

= 8.47 x 10~ 14 aA 2 e- |z|/A peV (14) 

with p = 19 g/cm 3 (gold or tungsten coating) and A given 
in pm. Taking these gravity-like forces into account, first 
order perturbation theory predicts a shift of the n-th en- 
ergy eigenvalue [22l |. 

AE n = ^ n \AV{z)\^ n ). (15) 

They differ from state to state in the range of interest. 
The sensitivity can been seen in Fig. [51 where the energy 
shift as a function of range A is plotted for a fixed a — 
10 12 . 

Newtonian gravity and hypothetical fifth forces 
evolve with different phase information in the non- 
oscillating region. We expect the following sensitivity for 
50 days of beam time at the PF2-UCN beam position at 
the ILL: With a count rate of 0.1 s _1 for neutrons in the 
ground state, we will have N = 430000 registered neu- 
trons. Due to the uncertainty principle A<i>A7V > 2ir, we 
estimate a minimal detectable phase shift of 9.6xl0~ 3 ra- 
dians. For an estimate of T = 130 ms interrogation time 
(flight path between the oscillators), the minimal resolv- 
able energy shift is 

AE = A$fi jT = 0.096/i/s = 4.8 x 10" 5 peV. (16) 

Together with Eq. (|14p. this corresponds to a sensitivity 
of a = 7.4 x 10 7 , which is about three orders of magni- 
tude better than existing limits right now. In principle, 
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FIG. 6: The relative frequency shift vs. range A. The strength 
of the modified Newtonian potential AV(z) of Eq.[l4]is set to 
a = 10 12 . 



ultracold neutrons can be stored and the time scale can 
be increased to T = 130 s. Therefore, a sensitivity to 
energy differences of AE — 4.8 x 10 -20 eV is feasible. 
This corresponds to an a = 7.6 x 10 4 . With new neu- 
tron sources, which are under development right now, the 
source strength density is expected to be increased by two 
orders of magnitude. The statistical sensitivity of the new 
method is therefore around 



or 



AE = 4.8 x l(T 21 eV 



a < 7.6 x 10 3 . 



(17) 



(18) 



This is orders of magnitude better than existing limits 
right now. 

Limits for hypothetical fifth forces can be easily in- 
terpreted as bounds of the strength of the matter cou- 
plings of axions. Axion interactions with a range within 
20 \im < A < 200 mm (corresponding to axion masses 
10 _6 eV < m a < 10~ 2 eV), the 'axion window', are still 
allowed by the otherwise stringent constraints posed by 
cosmological data (see e.g. [H, HH). The CP-violating 
spin-dependent part in presence of matter given by [57| is 



V(r) = hg p g s 



87rmc 



1 

a7 



i 



-r/X 



(19) 



related to the geometry of the macroscopic matter con- 
figuration. 

I nteg rating this potential over the geometry of region 
three [22j, we arrive at A = 5 /im at a limit of 



^ < 5.3 • 10- 23 
he 



(20) 



Here, a denotes the neutron spin and n is a unit vector 



for the dimensionless axion coupling strength. This is 
again eight orders of magnitude better than the only 
existing limit [22l . [23| in the axion window from the 
previous experiment with neutrons. 



IV. SUMMARY 



In conclusion, we discussed an application of Ram- 
sey's method of oscillating fields to the quantum bouncer. 
It will allow high precision spectroscopy of the energy 
eigenstates of a neutron bouncing on a flat vertical sur- 
face. Such Ramsey type interference measurements will 
improve the sensitivity for neutron's coupling to gravity, 
to hypothetical short ranged forces or the influence of the 
cosmological constant. A sensitivity of 21 orders of mag- 
nitude below the strength of electromagnetism is found, 
when the energy AE = 4.8xl0 -21 eV of Eq. (17) is com- 
pared with the Rydberg energy of 13.6 eV, which is the 
energy scale of electromagentically bound quantum sys- 
tems. Such an energy change corresponds to a strength 
a ~ 7.6 x 10 3 compared to gravity or to ~ 5.3 • 10~ 23 , 
the axion coupling strength, at a range A = 5 /an. 

The new method profits from small systematic effects 
in such systems, mainly due to the fact that in contrast 
to atoms, the electrical polarizability of neutrons is ex- 
tremely low. Neutrons are not disturbed by short range 
electric forces such as van der Waals or Casimir forces. To- 
gether with its neutrality, this provides the key to a sen- 
sitivity of several orders of magnitude below the strength 
of electromagnetism. 
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